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Suppose  a complex  mechanism,  e.g.,  a missile,  is  built 
from  a number  of  different  types  of  components,  where  the 
reliability  of  each  of  the  components  has  been  estimated  by 
means  of  separate  tests  on  each  of  the  components.  There 
exist:  many  procedures  for  combining  such  component  data  to 
determine  approximate  lower  confidence  limits  for  the 
reliability  of  the  system.  Several  such  procedures  are  the 
Maximum  Likelihood  (Ref.  12) , the  Madansky  (Ref.  5) , the 
Log-Gamma  (Ref.  13) , the  Easterling/Modified  Maximum 
Likelihood  (Ref.  2)  and  the  Mann  (Ref.  6)  methods. 

This  study  is  concerned  with  series  system.  Lower 
confidence  limit  accuracies  are  compared  by  means  of  computer 
simulation.  Two  methods  of  introducing  partial  component 
failures  were  used  in  the  simulation  to  modify  some  of  the 
procedures  which  cannot  be  used  when  all  the  components 
exhibit  no  failures.  The  first  method  was  proposed  by 
W.  M.  Woods  and  the  second  by  Lisowsky  (Ref.  4)  who  also 
developed  a successful  technique  for  computing  the  Lagrange 
multiplier  in  the  Madansky  procedure.  An  attempt  was  also 
made  to  improve  the  performance  of  the  Easterling  method 
using  randomization  techniques  developed  by  D.  R.  Barr  and 
T.  Jayachandran  (Ref.  1) . 

The  simulation  results  show  that  the  accuracy  of  the 
Maximum  Likelihood,  Madansky  and  Mann  procedures  are  comparable. 
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the  Log-Gamma  and  Easterling  procedures  yield  satisfactory 
results  when  sample  sizes  are  large  and  unequal,  and  the 
Randomized  Easterling  procedure  yields  better  results  than 
the  Easterling  procedure  in  all  cases. 


Finally,  it  is  worthwhile  to  note  that  the  Maximum 
Likelihood  procedure  is  simple  and  easy  to  implement; 
therefore  it  can  be  used  as  a rough  and  ready  method. 


II.  DESCRIPTION  OF  THE  METHODS 


Consider  a system  which  consists  of  k components  in 
logical  series;  the  components  may  be  either  continuously 
operating  or  of  the  cycle  type. 

Suppose  n^  copies  of  component  i are  put  on  test,  i = 1, 
2,...,k,  under  the  environmental  conditions  defined  in  its 
mission  profile,  and  let  each  operate  until  failure  or  the 
mission  time  is  reached,  whichever  occurs  first.  Denote  f^ 
as  the  number  of  components  of  type  i that  did  not  complete 
the  mission,  and  define 

pi  = 1 - qi  (2.1) 

where 

q±  . fi/ni 

The  following  methods  were  modified  in  Chapter  III, 
section  C,  since  some  of  the  procedures  cannot  be  used  when 
all  components  exhibit  no  failure,  i.e.,  the  Madansky,  the 
Easterling  and  the  Randomized  Easterling  procedures. 

A.  THE  MAXIMUM  LIKELIHOOD  (ML)  METHOD 

The  maximum  likelihood  estimate  for  system  reliability 
is 

* k „ 

R = H p.  (2.2) 

s i=l  1 

This  estimator  is  asymptotically  normal  in  distribution,  and 
its  variance  is  estimated  by  (Ref.  9) : 


(2.3) 


42  ■=  Rg 


lni-Xl) 

niXi 


where  x^  = n^-f^ 

The  ML  100(l-a)%  Lower  Confidence  Limit  (LCL)  for  system 
reliability  is  given  by 


Rs  ' L ( a ) = Rs 


z,  a 
1-a 


(2.4) 


where  z1_a  is  the  100 (1-a)  percent  point  of  the  standard 
normal  distribution. 


B.  THE  MADANSKY  (MD)  METHOD 

The  Madansky  method  is  based  on  the  well-known  results 
due  to  Wilks  (Ref.  11)  that  -2  lnp  is  distributed  asymptote 
ically  as  a chi-square  random  variable  with  one  degree  of 
freedom. 

P is  the  likelihood  ratio  test  statistic  given  by 
(Ref.  5) . 


where  BIN  (xi;pi)  = (XX)  P^i  (1_p  } ni_Xi 

The  numerator  is  maximized  under  the  additional  constraint 


the  denominator  is  an  unconstrained  maximization.  Values  of 


R included  in  the  two  sided  confidence  interval  (the 
s 

confidence  set)  are  given  by 


s (R  ) = [R  : - 2 lnp  < xf  nl  (2.6) 

2 

where  x 1 is  the  upper  percent  point  of  the  chi-square 
distribution  with  1 degree  of  freedom. 

It  is  easy  to  see  that  the  logarithm  of  the  denominator 
of  p is 


k x.  k n.  k k 

I In  BIN(x.;— =■)  = E In  ( :)  + I x.lnx.  - E x.ln  n. 
i=l  1 ni  i=i  xi  i=l  1 1 i=l  1 1 


k k 

+ £ (n.-x.)  In  (n.-x.)  - £ (n.-x.)  In  n. 

i=l  1 1 11  i=l  1 1 1 


To  determine  the  logarithm  of  the  numerator  of  P , let  us 

first  maximize  the  Lagrangian 

k k 

£ In  BIN  (x . ; p.)  - A {In  Ft  p.  - In  R 1 
i=l  1 1 i=l  1 s 


where  X is  a Lagrange  multiplier.  The  maximizing  set  of 
p^'s  is  given  by 


where,  since  0 <_  <_  1 for  all  i,  then  X <_  min  x^. 

Hence,  as  a function  of  A, 


k n. 


E In  BIN  (x.  ; p. ) ■ l (v  ) 
i=l  1 1 i-1  Xi 


k k 

+ Z x.  ln(x.-X)  - Z x,  ln(n.-X) 


i-i'1  1 ' i-i 1 


fw 

+ E (n. -x. ) In (n. -x. ) - E (n.-x. ) ln(n.-X) 
i=l  1 1 11  i=l  1 1 


and 


lnp 


E x.  In  (1  - - Z n.  In  (1  - -£-)  , 

i=l  1 xi  i=l  x i 


(2.7) 


where  X satisfies  the  constraint  equation 


k (xi-X) 


n , ---T-  = r 

i-i  s 


(2.8) 


The  set  of  X such  that  -2  In  p < x n will  be  an  interval 

a , x 

* *.  * * * . , 

(Xr  X2]f  where  X1  < 0 < X2  and  the  X satisfy 

2 


E x.  ln(l  - - 2 n.  In  (1  - ~)  r 


- -L)  = - 


i=l 


xi  i=l  1 


n . 

i 


(2.9) 


The  MD  100  (1-a) % LCL  is  given  by 

* 


(2.10) 


Rs'L(a)  = Rs(X2} 

In  determining  lower  confidence  limits  the  convention  adopted 
by  Myhre  and  Saunders  (Ref.  12)  will  be  used  where  the  value 
of  a in  equation  2.9  is 

o = 2 (1-y)  , 


and  y is  the  confidence  coefficient 


C.  THE  LOG-GAMMA  (LG)  METHOD 

In  the  LG  procedure  the  method  of  moments  is  used  to  fit 
the  random  variable  -InR  with  the  two  parameter  gamma 
distribution  (Ref.  13) . The  gamma  is  then  transformed  into 
a chi-square  distribution,  about  which  probability  statements 
are  made  and  a lower  confidence  limit  obtained.  That  is 
define 

k 

S * -InR  = - Z In  (1-q . ) (2.11) 

s i-1  1 

Expanding  the  natural  logarithm  in  an  infinite  series 

S « - Z [(-q.)  - 3s  (-q  • ) 2 + i •••] 
i=l  1 1 


Mil 


1 


2n.  - 3 

where  a.  =*  w-7-i , 

1 2(ni  - : 


(2.14) 


b.  = — 

1 ni  - 


That  Ti  is  unbiased  is  important,  because 


(2.15) 


- S a Z T. 


(2.16) 


is  used  as  an  estimator  for  S,  thereby  accumulating  any  bias 
present  in  the  . An  approximate  value  for  the  variance  of 
S (Ref.  7)  is 


k k T 

Var  (S)  = Z Var  (T.  ) = z -i- 
i=l  1 i=l  ni 


(2.17) 


* 

Next,  fit  S with  a gamma  distribution.  The  probability 
distribution  of  S is  then  given  by  the  density  function 


f (x;  r,  0)  = xr_1  exp  (-x/9) 

r (r)  9r 


(2.18) 


x>0,  r>0,  9>0 


It  follows  that 


E (S)  = re 
var  (S)  = re2 

Since  S is  unbiased 
E (S)  « S 

k 

= Z T. 
i=l  1 


(2.19) 

(2.20) 


(2.21) 


Solving  equations  (2.17),  (2.19),  2.20)  and  (2.21) 
simultaneously  gives  the  shape  parameter 


MA  * 


and  the  scale  parameter 


fk  TI1  A 

li=i  w i=1 


Thus,  r can  be  estimated  by 

2 

Is  * I 

r 


“k  A“12  / k T. 
2 T . / S — 

i-1  / i-1  ni 


(2.22) 


(2.23) 


(2.24) 


Since  S is  distributed  gamma  (r,0),  2S/9  is  distributed 
2 2 

x2r*  Then  X1_a  2r  is  that  number  such  that 


1 - ° ■ p t;*2r  i *l-„,2r 
And  since  S is  unbiased 
E (S ) = S 
= r9 

= _lnRs 

Equation  2.25  becomes 


(2.25) 


(2.26) 


1 - « " P <2rS  .>  -lnRs  xl-a,2rJ 

-2rS. 

,2r . 


- a = P F"exp  l-^2ri 

L 


^Rs 


(2.27) 


2 


Define  as 


(f±  + 1) 

o 


+ % bi 

o 


(f±  + l)2 

o 


This  means  that  one  more  failure  is  added  to  that  component 

I 

with  largest  sample  size  to  obtain 

A 

3.  Define  T*  by 


T*  = k (T!  + T.  ) 

i l i ' 

o o o 

4.  Substitute  T7  for  in  S to  obtain  S* 

o a o 

The  resulting  S*is  the  continuity  corrected 
value  of  S. 

5.  r is  corrected  to  obtain  r*  by  substituting  T£ 

A A 

for  T.  in  the  definition  of  r. 

10  A A 

6 . With  these  definitions  of  S*  and  r*  the  new 

A 

100  (l-a)%  lower  confidence  limit  R*  r , « for 

s ,L  (a) 

R becomes 
s 


Hs ,L(a) 


exp 


- [2r*]  S* 


Xl-a,{2r*] 


(2.31) 
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D.  THE  EASTERLING /MODIFIED  MAXIMUM  LIKELIHOOD  (MML)  METHOD 

In  this  method  the  ML  estimate  R is  treated  as  the 

s 

usual  binomial  estimate  based  on- "n",  called  the  pseudo  sample 
size,  is  unknown  and  is  estimated  from 


a2 


Rs(1  - 


V 


n 


(2.32) 


* 2 

where  a is  given  by  2.3.  Thus  by  equating  the  estimated 

variance  of  R under  maximum  likelihood  theory  to  what  it 
s 

would  be  under  binomial  theory,  we  can  solve  this  equation 
for  n.  Then  the  component  test  results  can  be  regarded  as 

A Ai  A A 

being  equivalent  to  system  results  of  n tests  with  x = R .n 
successes . 

In  binomial  sampling  with  x successes  in  n trials,  a 
lower  100(1  - a) % confidence  limit  on  the  reliability  is 
given  by  the  solution  for  Rg  in 


where 


L 


Tu 


(R  , x,  n - x + 1)  = « 

S 


(2.33) 


the  incomplete  beta  function  with 
x : the  first  parameter,  and 
n - x + 1 : the  second  parameter, 
n and  x are  unlikely  to  be  integers  and  the  calculations  of 
lower  limits  in  the  comparative  examples  (chapter  IV)  use 
Easterling's  MML  method  in  which  n and  x were  rounded  up  to 
the  next  integers. 
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E.  THE  RANDOMIZED  EASTERLING  (RE)  METHOD 

Barr  and  Jayachandran  in  Ref.  1 develop  a randomization 
technique  for  improving  the  lower  confidence  limit  on  the 
reliability  of  a system  with  a discrete  distribution. 

Basically  the  method  is  a simple  one,  for  example  an  exact 
100(1  - a)%  LCL  for  the  parameter  of  the  binomial  distribution 
can  be  obtained  as  follows : 

1.  Let  z=x+y,  where  y represents  an  observed  value  of 
Y,  and  Y is  uniformly  distributed  between  zero  and  one 
inclusive;  x is  the  number  of  successes  in  n trials. 

2.  The  solution  for  Rs  in 

_Z~(RS/  z,  n - z + 1)  « a (2.34) 

is  the  exact  100(l-a)%  LCL. 

Since  the  Easterling  method  leads  to  one  component 
system,  equation  2.34  is  readily  applicable  for  computing  an 
exact  100(l-a)%  LCL  based  on  the  Easterling  method  by 
replacing  z with  z in  equation  2.35,  where 
z = x + y 

and  n with  n,  x and  n as  defined  in  equation  2.33. 

F.  THE  MANN  (MN)  METHOD 

Mann  et  al.  in  Ref.  6 shows  that  the  100(l-a)%  LCL  on 

R can  be  obtained  by  fitting  the  posterior  distribution  of 
s 

-lnRg  with  a noncentral  chi-square  distribution.  The 
corresponding  central  chi-square  variate  with  non-integer 
degrees  of  freedom  is  transformed  to  normality  yielding: 


m 


(2.35) 


where  m is  the  mean  and  v the  variance  of  the  posterior 
distribution  of  -lriRs. 


m 


0.5 (1+1/a) 

3 


1-RS 

0.5 (Rs+i) 


(2.36) 


with 


n 


(2.37) 


(2.38) 


n^j  = min (n1/ . . . ,nk)  , (2.39) 

n*  = n{1)  [1-0.5(1-RS)2]  [1-0. 5(  1-Rs)  ] , (2.40) 

and  z1_a  “ 100 (1-a) -th  percentile  of  the  standard  normal 
distribution. 

All  zero  - failure  components  are  ignored  in  calculating 
the  value  of  a in  equation  2.38  except  for  any  single  zero  - 
failure  component  with  its  sample  size  equal  to  n^j.  It, 
too,  is  ignored,  however,  if  at  least  one  other  component 
that  exhibits  failure  has  sample  size. equal  to  n^. 
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III.  THE  SIMULATION 


A.  GENERAL  METHODOLOGY 

A computer  simulation  is  used  in  this  study  as  an 
analytical  tool  for  evaluating  a proposed  LCL  procedure;  the 
method  is  as  follows: 

1.  Suppose  it  is  desired  to  evaluate  a proposed 
100  (1-a) % LCL  procedure,  denoted  R . , for  system  reli- 

s , li  va  / 

ability  Rg.  Then  the  assertion  is 


t _ , . < R ] > 1-a 
s,L(a)  — s — 


(3.1) 


Equality  should  hold  if  R . , , is  a continuous  random 
variable . 

2.  Rg  = f (p1,p2, . . . ,pk) , where  pi  is  the  true  relia- 
bility of  the  component.  For  an  independent  series  system, 

k 

Rs  - n Pi* 
s i=l  1 

3.  Assign  values  to  the  parameters  a,  k,  n^  and 
p^,  i=l,2,...,k.  Perform  n^  Bernoulli  trials  for  the  ith 
component  to  get  the  number  of  successes/failures  for  each 
component  and  then  compute  the  resultant  Rg  L(a)* 

4.  Generate  the  approximate  distribution  of  Rg  L(a) 
by  repating  step  3 500  times. 

5.  Order  the  R„  _ , . realizations  to  get 

s ,L (a) 

, A A 

Rs  /L  (a)  s,L(a)  (50o)'’ 


6.  Find  the  500(l-a)th  order  statistic  of  R T,  , 

— s ,i,  la ) 

and  denote  it  A^_a>  Thus  A^_a  is  the  (1-g) th  percentile 

A 

of  the  distribution  of  the  LCL  random  variable  R T , N . 

s ,L (a) 

7.  By  repeating  steps  3 through  6 for  various  sets  of 
k/n^^r . . . ,n^,  Pi'P2'  * * * 'Pjc  comparing  the  resultant 
A1_o  to  Rg,  the  overall  performance  of  the  LCL  procedure  can 
be  evaluated. 

The  actual  confidence  level  given  by  R T . s can  be 

s ,L  (a; 

obtained  by  finding  the  order  statistics  of  the  generated 
distribution  which  matches  (or  is  closest  to  matching)  Rg . 

it 

If  the  index  of  this  order  statistic  is  denoted  i , then 


x 100%  = Actual  level  of  confidence 


(3.2) 


of  R 


s,L(a) 


Equivalently,  if  < Rg , the  procedure  is  a conservative 

one,  and  vice-versa. 


B.  ACCURACY  CRITERIA 

There  are  three  characteristics  of  the  distribution  of 

R_  _ , » for  determining  the  accuracy  of  LCL  procedure, 
s , Xj  \ a ; 

They  are: 

1.  The  mean 

2.  The  variance,  and 


3.  A 


1-ct* 


The  variance  should  be  small  and  the  actual  values  should 
still  be  within  the  ball  park  when  the  LCL  is  applied. 


p ' 


If  in  fact,  T / v is  an  exact  100(l-a)%  LCL  procedure 

S f Li  \ Cl  ) 


for  Rg , that  is 


(3.3) 


P(R  , . < RJ  = l-o 

then  A.  should  be  close  to  R„  regardless  of  the  set  of 
l— a s 

parameter  values  used  and  for  each  value  of  a. 

Thus  the  quantity* 


Al-a  "Rs|  (3*4) 

is  a measure  of  the  accuracy  of  the  procedure. 

C.  SIMULATION  ALGORITHM 

Step  1.  Given  a set  of  parameters  k,  n^  and  p^, 
i * l,2,...,k,  generate  binomial  data  as  follows: 

For  the  ith  component,  i = l,2,...,k,  draw  a uniformly 
distributed  random  number  u from  the  interval  [0,1]  and 
compute  the  number  of  successes  x^ 


si  =< 


if  u < Pj_ 
,0  if  u > Pi 


n. 

i 

x • = 2 s . , x = 1,2, ...,k. 

i ^ i 


(3.5) 


All  the  procedures  described  in  chapter  II  simply 
ignore  the  zero-failure  component (s)  with  the  result  either 
the  computed  LCL  tends  to  be  very  high  or  the  procedures 
cannot  be  applied. 

2 

For  example  if  for  all  i,  then  a in  equation  2.3 

is  zero,  and  hence  the  ML's  LCL  in  equation  2.4  is  equal  to 
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one  and  the  Easterling  -nd  Randomized  Easterling  procedures 

cannot  be  applied  because  n,  the  pseudo  sample  size  in 
equation  2.32,  is  not  defined.  Likewise  the  Madansky  proce- 
dure cannot  be  applied  (equation  2.9  is  not  defined) . In 
this  case  the  Log-Gamma's  LCL  and  the  Mann's  LCL  are  close 
to  one.  To  cope  with  this  problem  two  methods  of  introducing 
partial  component  failures  will  be  presented  in  this  study 
for  those  methods  which  accommodate  zero  failures. 

FIRST  METHOD.  If  x^sn^  for  all  i,  pick  the  component 
with  the  largest  sample  size  or  pick  the  first  component  if 
all  sample  sizes  are  equal  and  introduce  a half  failure  to 
that  component;  this  method  was  proposed  by  W.  M.  Woods. 

SECOND  METHOD.  Whenever  fi=0 , set  fi=4/(nik),  this 
method  was  proposed  by  Lisowsky  (Ref.  4) . 

Step  2.  For  each  a (asO.l  and  a^0.2)  compute  the  LCL 
for  all  procedures  as  described  in  chapter  II . 

Step  3.  For  each  set  of  data  repeat  step  1 and  step  2 
500  times,  order  the  LCL ' s , and  then  compute  A^_a , mean, 
standard  deviation  of  the  LCL  distribution,  the  actual  level 
of  confidence  in  equation  3.2,  and  the  average  failure  per 
replica. 

Now  determine  each  procedure's  accuracy  by  comparing 
the  Ai_a's  with  the  true  system  reliability  Rg,  and  relative 
dispersion  by  comparing  their  mean  and  standard  deviation. 

For  detailed  computation  algorithm  see  Appendix,  The 
Computer  Program. 
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IV.  SIMULATION  RESULTS  AND  CONCLUSIONS 


A.  GENERAL  DESCRIPTION 

The  accuracy  of  the  procedures  are  compared  for  a variety 
of  sets  of  parameter  values  (k,  pi,  n^,  i=l,2, . . . ,k) . These 
different  sets  of  parameters  are  called  cases  and  are 
numbered.  For  each  case  the  two  methods  of  introducing 
partial  component  failures  were  applied  at  two  different 
confidence  levels:  90%  and  80%  CL. 

The  simulation  results  were  tabulated  in  Table  II  at 
the  end  of  this  chapter;  for  each  case  the  results  were 
listed  in  one  table,  for  example  Table  II. 1 for  case  number 
1,  Table  II. 2 for  case  number  2,  etc.,  up  to  Table  11.36. 

The  letter  suffix  attached  to  the  case  number  denotes  the 
method  of  introducing  partial  component  failures.  In 
particular  case  la  means  the  FIRST  METHOD  of  introducing 
partial  component  failures  as  described  on  page  24  was  applied 
to  case  1.  Similarly  case  lb  means  the  SECOND  METHOD  of 
introducing  partial  component  failures  as  described  on  page  24 
was  applied  to  case  1. 

By  varying  n^  while  k and  p^  are  held  constant  in  each 
case,  the  sensitivity  of  the  procedures  can  be  examined. 

B,  THE  EFFECT  OF  SAMPLE  SIZE 

Observing  the  simulation  results  in  Table  II  starting 
from  sample  size  15  through  50  on  each  case,  these  results 
indicate  that  increasing  sample  size  (continued  on  page  27  ) 


TABLE  1 


k 
n . 


m 

s 

Fbar 


Nfail 

Ncorr 

CL 

ACL 

A. 

1-a 

ML 

MD 

LG 

MML 

RE 

MN 


NOMENCLATURE  FOR  SIMULATION  RESULTS 

Number  of  series  connected  components 

Sample  size  of  ith  component 

True  reliability  of  ith  component 

True  independent  series  system  reliability, 

k 

R = n P4 
S i=l  1 

Sample  mean  of  R T / \ 
r s ,L  (a  ) 

Sample  standard  deviation  of  R _ , . 
c s ,L  (a ) 

Average  failure  per  replica 
500  k ni 

(1/500)  E Z E f. 

1 i=l  1 1 

Total  number  of  FIRST  METHOD  of  introducing 
partial  component  failures  applied 

Total  number  of  SECOND  METHOD  of  introducing 
partial  component  failures  applied 

100 (1-a) % confidence  level 

Actual  confidence  level,  computed  using  equation 
3.2,  page  22,  the  result  is  rounded  to  the 
smallest  integer 


The  500 (1-a) th  order  statistic  of  R 
Maximum  likelihood  method 
Madansky  method 
Log-Gamma  method 
Easterling  method 
Randomized  Easterling  method 
Mann  method 


s , L ( a ) 
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causes  the  A 


l_a's  to  converge  to  Rg,  thus  decreasing  the  value 
of  |A]__a"Rs|  of  all  procedures.  The  results  also  indicate 
that  some  of  the  procedures  need  larger  sample  sizes  to 
converge  satisfactorily. 

With  medium  sample  size  (30)  and  large  sample  size  (50) 

the  performances  of  ML,  MD  and  MN  procedures  are  comparable, 

although  at  sample  size  30  in  some  cases  MD  and  MN  procedures 

perform  better  than  ML  procedure  as  shown  in  case  5,  14  and 

23.  But  with  sample  size  15  and  with  R > 0.900  the  ML 

s 

procedure  seems  to  perform  better  than  MD  and  MN  procedures 
as  shown  in  cases  1,  16  and  19. 

The  LG,  MML  and  RE  procedures  converge  rather  slowly, 
but  RE  tends  to  converge  faster  than  MML  in  all  cases;  this 
is  due  to  the  randomization  effect  as  described  in  chapter  II 
section  E page  .19.  On  the  other  hand  with  unequal  sample 
sizes  (mix  of  large,  medium  and  small  sample  sizes)  these 
procedures  tend  to  converge  satisfactorily  as  shown  in  case 
34,  35  and  36,  while  the  MN  procedure  always  yields  unsatis- 
factory results,  and  again  ML  and  MD  procedures  are  still 
comparable  under  this  condition.  Note  that  in  these  cases 
only  the  SECOND  METHOD  of  introducing  partial  component 
failures  can  be  applied  in  the  MD  procedure  as  was  pointed 
out  by  Lisowsky  (Ref.  4) . 


C.  THE  EFFECT  OF  THE  METHOD  OF  INTRODUCING  PARTIAL 

COMPONENT  FAILURES 

When  the  cases  under  study  have  large  sample  sizes  if 
it  can  be  expected  that  all  components  of  these  cases  exhibit 
failure (s)  or  only  a few  components  exhibit  no  failure,  then 
the  two  methods  practically  have  no  effect  on  the  computed 
LCL’s  as  demonstrated  by  simulation  results  using  sample 
sizes  of  50  in  Table  II.  These  facts  were  also  shown  in  cases 
where  p^  < 0.960  for  all  i with  sample  size  30,  i.e.,  case  8, 

11  and  14,  Table  II.  The  explanation  is  as  follows.  Take 
as  an  example  case  8a  where  Nfail  = 2;  this  means  that  the 
FIRST  METHOD  of  introducing  partial  component  failures  was 
applied  twice.  Its  contribution  to  the  average  failure  per 
replica  (Fbar  = 5.6)  was  very  small:  (2)  (0.5) (1/500)  = 

0.002.  In  case  8b  Ncorr  = 768,  the  SECOND  METHOD  was  applied 
768  times  and  its  contribution  of  failure  to  Fbar  (5.7)  was 
{( 768x4) / (5x30) } (1/506)  = 0.04.  Both  of  these  contributions 

were  small  compared  to  Fbar  so  that  the  computed  A1_a's  in 
case  8a  and  8b  were  approximately  the  same.  Note  that  the 
value  of  Fbar  is  rounded  to  the  first  decimal  in  Table  II. 

On  the  other  hand  if  n^  is  small  and  p^  is  high  for  all 

i as  in  cases  1,  16  and  19,  then  different  values  of  A. 

1-a 

can  be  expected  from  the  two  methods.  For  example  case  la, 

Nfail  = 251  and  the  contribution  of  failure  to  Fbar  (1.0) 

I 

was:  (251) (0.5) (1/500)  = 0.25.  In  case  lb  Ncorr  = 2171, 

its  contribution  to  Fbar  (0.9)  was  { (2171x4)  / (5x15)}  (1/500)  = 

0.23.  Both  of  these  contributions  were  quite  high  compared 


with  Fbar  and  therefore  the  computed  A^_o's  were  different, 

except  for  MML  and  RE  procedures  which  were  due  to  the  fact 

that  n and  x in  chapter  II  section  D and  E on  pages  18  and  19 

rounded  up  to  the  next  integer.  In  these  cases  the  ML,  MD 

and  MN  procedures  yield  higher  and  the  LG  procedure  yields 

lower  A,  's  with  the  SECOND  METHOD,  but  this  was  not  so  for 
l-o 

the  other  cases. 

A more  detailed  quantitative  analysis  is  needed  to  deter- 
mine the  range  of  values  of  k,  n^  and  p^  for  which  the  FIRST 
METHOD  is  better  than  the  SECOND  METHOD  for  each  procedure 
and  vice  versa.  This  is  beyond  the  scope  of  this  study. 


D.  CONCLUSIONS  AND  REMARKS 

The  overall  performance  of  ML,  MD  and  MN  procedures  are 
comparable  although  with  medium  sample  size  (thirty)  MD  and 
MN  procedures  perform  better  than  the  ML  procedure,  but  with 
small  sample  size  (fifteen)  the  ML  procedure  seems  to  perform 
better  than  MD  and  MN  procedures . 

The  LG,  MML  and  RE  procedures  tend  to  yield  satisfactory 
results  using  unequal  sample  sizes  (mix  of  large,  medium  and 
small  sample  sizes) , while  the  MN  procedure  always  yields 
unsatisfactory  results.  The  ML  and  MD  procedures  are  still 
comparable  under  these  conditions. 

The  randomization  technique  was  successful  since  RE 
procedure  yielded  better  results  than  MML  procedure. 
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The  ML  procedure  is  simple  and  easy  to  implement,  there- 
fore this  procedure  can  also  be  used  as  a rough  and  ready 
method.  The  MD  procedure  requires  a computer  to  work  with 
while  the  MN  procedure  can  be  solved  with  a handheld  calcula 
tor  although  it  is  a rather  complicated  computation. 


TABLE  II.  1 


TABULATED  SIMULATION  RESULTS 


Case 

la: 

k=5,  p. 

=0.990, 

n.*15, 

i=l , 2 , . 

• • / 5 • 

u 

o 

• 

951 

FIRST  METHOD:  Fbar  = 1.0,  Nfail  = 251 

********90%  CL**********  *********90% 

Cl********* 

A.90 

m 

s 

ACL 

A.80 

m 

s 

ACL 

ML 

.907 

.861 

.060 

50 

.928 

.887 

.054 

50 

MD 

.872 

.830 

.056 

50 

.912 

.874 

.052 

50 

LG 

.708 

.660 

.053 

50 

.818 

.776 

.049 

50 

MML 

.778 

.734 

.053 

50 

.824 

.782 

.051 

50 

RE 

.847 

.773 

.062 

100 

.873 

.818 

.059 

100 

MN 

.810 

.777 

.042 

50 

.855 

.823 

.041 

50 

Case  lb:  same  as  la 

SECOND  METHOD:  Fbar  = 0 

ML  .939  .869  .079 

.9,  Ncorr  = 2171 

50  .954  .893 

.070 

50 

MD 

.901 

.837 

.073 

50 

.938 

.879 

.067 

50 

LG 

.652 

.627 

.020 

82 

.797 

.766 

.039 

50 

MML 

.778 

.728 

.056 

50 

.824 

.777 

.054 

50 

RE 

.847 

.768 

.065 

100 

.873 

.814 

.061 

100 

MN 

.885 

.829 

.065 

50 

.914 

.862 

.061 

50 

TABLE  II. 2 (Continued) 


f 


r 

» 


Case  2a:  Same  as  la,  except  n^SO,  i=l,2,...,5 
FIRST  METHOD:  Fbar  = 1.6,  Nfail  = 121 

*******4rgQ%  0^**********  ****r****gQ£  0^********* 


A.90 

m 

s 

ACL 

A.80 

m 

s 

ACL 

ML 

.953 

.899 

.050 

76 

.964 

.915 

.045 

45 

MD 

.934 

.882 

.048 

76 

.956 

.908 

.044 

76 

LG 

.842 

.791 

.036 

76 

.905 

.857 

.038 

76 

MML 

.880 

.824 

.045 

76 

.906 

.854 

.042 

76 

RE 

.905 

.845 

.048 

100 

.913 

.873 

.045 

100 

MN 

.901 

.861 

.036 

76 

.925 

.887 

.035 

76 

Case  2b:  Same  as  2a 

SECOND  METHOD:  Fbar  = 1.6,  Ncorr  = 1844 


ML 

.980 

.902 

.058 

76 

.985 

.918 

.051 

45 

MD 

.959 

.885 

.055 

76 

.977 

.911 

.050 

51 

LG 

.819 

.797 

.029 

76 

.910 

.862 

.040 

76 

MML 

.880 

.823 

.044 

76 

.906 

.853 

.042 

76 

RE 

.905 

.843 

.048 

100 

.913 

.872 

.045 

100 

MN 

.948 

.884 

.050 

76 

.962 

.904 

.046 

76 

Case  3a:  Same  as  la,  except  n^=50,  i=l,2,...,5 
FIRST  METHOD:  Fbar  = 2.5,  Nfail  =49 

********90%  cl**********  ********3Q%  cl********** 


A.90 

m 

s 

ACL 

o 

00 

• 

< 

m 

s 

ACL 

ML 

.955 

.915 

.040 

68 

.963 

.928 

.037 

68 

MD 

.944 

.905 

.039 

86 

.959 

.924 

.036 

68 

LG 

.893 

. 860 

.032 

90 

.930 

.896 

.033 

90 

MML 

.901 

.868 

.036 

90 

.920 

.889 

.034 

90 

RE 

.926 

.880 

.038 

98 

.931 

.900 

.035 

93 

MN 

.925 

.897 

.031 

90 

.941 

.914 

.031 

90 

Case  3b:  Same  as  3a 

SECOND  METHOD:  Fbar  = 2.5,  Ncorr  = 1529 


ML 

.953 

.916 

.043 

68 

.962 

.928 

.039 

68 

MD 

.942 

.906 

.042 

72 

.957 

.924 

.038 

68 

LG 

.889 

.860 

.030 

90 

.927 

.897 

.034 

90 

MML 

.899 

.868 

.036 

90 

.918 

.888 

.034 

90 

RE 

.925 

.880 

.037 

98 

.929 

.900 

.035 

93 

MN 

.941 

.908 

.038 

73 

.953 

.921 

.036 

68 

TABLE  II. 4 (Continued) 


Case  4a:  k=5, 

FIRST  METHOD: 

p^*0.977, 

Fbar  = 1. 

n^=15,  i=l,2,. 

8,  Nfail  = 84 

• • f 5 • 

R =0 . 890 
s 

******** 

90% 

££********** 

********30% 

cl********** 

A.90 

m 

s 

ACL 

> 

• 

00 

o 

m 

s 

ACL 

ML 

.907 

786 

.096 

83 

.879 

.819 

.087 

53 

MD 

.872 

760 

.089 

83 

.866 

.808 

.084 

83 

LG 

.708 

614 

.057 

83 

.780 

.721 

.066 

83 

MML 

.778 

671 

.076 

83 

.765 

.720 

.074 

83 

RE 

.818 

706 

.083 

100 

.825 

.754 

.080 

97 

MN 

.810 

725 

.069 

83 

.813 

.771 

.069 

83 

Case  4b:  Same  as  4a 

SECOND  METHOD:  Fbar  - 1.9,  Ncorr  = 1754 


ML 

.939 

.780 

.103 

83 

.861 

.813 

.093 

65 

MD 

.901 

.754 

.095 

83 

.849 

.802 

.oyo 

53 

LG 

.652 

.609 

.044 

81 

.759 

.717 

.064 

83 

MML 

.778 

.665 

.074 

83 

.751 

.716 

.072 

83 

RE 

.818 

.701 

.081 

100 

.815 

.750 

.079 

97 

MN 

.885 

.754 

.089 

83 

.837 

.791 

.084 

83 

TABLE  II. 5 (Continued) 


Bn. 


Case  5a:  Same  as  4a , except  n^=30 , i=l,2,...,5 
FIRST  METHOD  : Fbar  = 3.5,  Nfail  = 10 

********90%  cl**********  ********30%  cl********** 


A.90 

m 

s 

ACL 

> 

• 

00 

o 

m 

s 

ACL 

ML 

.925 

.818 

.068 

75 

.897 

.842 

.063 

70 

MD 

.907 

.805 

.065 

87 

.890 

.836 

.062 

85 

LG 

.795 

.741 

.046 

98 

.848 

.796 

.052 

95 

MML 

.842 

.754 

.058 

98 

.834 

.787 

.056 

98 

BE 

.853 

.773 

.061 

98 

.856 

.805 

.1)58 

92 

MN 

.877 

.799 

.055 

98 

.877 

.8  26 

.054 

75 

Case  5B  : Same  as  5a 

SECOND  METHOD:  Fbar  =3.6,  Ncorr  = 1217 


ML 

.919 

.816 

.068 

76 

.893 

.840 

.063 

70 

MD 

.902 

.803 

.065 

87 

.887 

.834 

.062 

83 

LG 

.818 

.744 

.050 

98 

.845 

.797 

.055 

89 

MML 

.837 

.753 

.057 

98 

.834 

.787 

.055 

98 

BE 

.849 

.772 

.060 

98 

.856 

.804 

.057 

93 

MN 

.902 

.808 

.062 

87 

.883 

.832 

.059 

76 

7* 


TABLE  II. 6 (Continued) 


Case  6a:  Same  as  4a,  except  n^=50,  i=l,2,...,5 
FIRST  METHOD:  Fbar  = 5.7,  Nfail  = 1 


***** 

***90% 

QJj **  ******** 

**********gQ^ 

Cl  ******** 

A.90 

m 

s 

ACL 

o 

00 

• 

< 

m 

s 

ACL 

ML 

.899 

.838 

.052 

82 

.891 

.856 

.049 

71 

MD 

.889 

.830 

.051 

91 

.887 

.852 

.048 

81 

LG 

.851 

.796 

.045 

98 

.864 

.831 

.045 

89 

MML 

.855 

.789 

.047 

98 

.855 

.823 

.045 

92 

RE 

.872 

.809 

.048 

95 

.874 

.833 

.045 

89 

MN 

.890 

.832 

.047 

87 

.886 

.850 

.045 

82 

Case 

6b: 

Same  as  6a 

SECOND  METHOD:  Fbar  « 5. 

7,  Ncorr  = 779 

ML 

.898 

.837 

.052 

82 

.891 

.855 

.049 

79 

MD 

.889 

.829 

.050 

90 

.887 

.852 

.048 

79 

LG 

.851 

.796 

.045 

98 

.864 

.830 

.045 

90 

MML 

.855 

.798 

.047 

98 

.855 

.822 

.045 

92 

RE 

.872 

.809 

.047 

95 

.874 

.833 

.045 

89 

MN 

.892 

.835 

.049 

82 

.886 

.852 

.047 

81 

36 


TABLE  II. 7 (Continued) 


rw 


Case  7a:  p. =0.961,  n.=15,  i=l,2,...,5  R =0.820 

X ii  S 

FIRST  METHOD:  Fbar  = 3.0,  Nfail  = 25 

*********90%  CL*********  ********80%  CL********** 


A.90 

m 

s 

ACL 

o 

00 

• 

< 

m 

s 

ACL 

ML 

.851 

.700 

.115 

79 

.879 

.740 

.107 

63 

MD 

.820 

.680 

.107 

79 

.866 

.731 

.104 

78 

LG 

.631 

.570 

.057 

95 

.742 

.662 

.082 

95 

MML 

.716 

.605 

.088 

95 

.765 

.655 

.087 

95 

RE 

.766 

.636 

.094 

97 

.772 

.686 

.092 

91 

MN 

.765 

.660 

.088 

95 

.813 

.706 

.089 

79 

Case  7b:  Same  as  7a 

SECOND  METHOD:  Fbar  = 3.1,  Ncorr  = 1381 


ML 

.831 

.692 

.115 

79 

.861 

.732 

.107 

76 

MD 

.802 

.673 

.106 

82 

.849 

.724 

.103 

79 

LG 

.652 

.567 

.065 

91 

.755 

.657 

.082 

95 

MML 

.700 

.601 

.085 

95 

.751 

.652 

.084 

95 

RE 

.753 

.633 

.091 

97 

.762 

.683 

.090 

94 

MN 

.801 

.677 

.103 

84 

.837 

.717 

.099 

79 

TABLE  II. 8 (Continued) 


Case  8a:  Same  as  7a,  except , n^=30 , for  all  i,  i=l,  2,..., 5 
FIRST  METHOD:  Fbar  = 5.6,  Nfa’il  = 2 

*********90%  CL**********  ********80%  CL********** 


A.90 

m 

s 

ACL 

A.80 

m 

s 

ACL 

ML 

.835 

.742 

.078 

83 

.823 

.771 

.074 

69 

MD 

.822 

.732 

.075 

84 

.818 

.766 

.073 

83 

LG 

.764 

.683 

.061 

99 

.784 

.734 

.065 

89 

MML 

.772 

.689 

.067 

91 

.771 

.724 

.066 

91 

RE 

.801 

.706 

.070 

93 

.794 

.741 

.068 

86 

MN 

.821 

.734 

.070 

87 

.815 

.762 

.068 

83 

Case  8b:  Same  as  8a 

SECOND  METHOD:  Fbar  =5.7,  Ncorr  = 768 


ML 

.833 

.741 

.078 

83 

.822 

.769 

.073 

71 

MD 

.820 

.731 

.075 

88 

.817 

.765 

.072 

81 

LG 

.762 

.684 

.064 

99 

.783 

.734 

.066 

89 

MML 

.772 

.689 

.067 

91 

.771 

.724 

.065 

91 

RE 

.800 

.706 

.069 

94 

.794 

.740 

.067 

86 

MN 

.825 

.738 

.073 

83 

.816 

.765 

.070 

82 
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TABLE  II. 9 (Continued) 


Case  9a:  Same  as  7a  except  n^=50  for  all  i,  i*l,2,...,5. 
FIRST  METHOD:  Fbar  =9.6  Nfail  = 0 

*********93%  0^*********  ********gg%  cl********** 


A.90 

m 

s 

ACL 

> 

• 

CO 

c 

m 

s 

ACL 

ML 

.829 

.757 

.060 

84 

.827 

.779 

.057 

76 

MD 

.821 

.751 

.058 

85 

.824 

.111 

.057 

76 

LG 

.790 

.725 

.054 

95 

.804 

.759 

.054 

84 

MML 

.792 

.726 

.054 

95 

.796 

.752 

.053 

87 

RE 

.808 

.736 

.055 

92 

.803 

.762 

.054 

86 

MN 

. 826 

.757 

.057 

85 

.824 

.778 

.055 

76 

Case  9b:  Same  as  9a 

second  METHOD:  Fbar  = 9.6,  Ncorr  = 376 


ML 

.828 

.757 

.060 

84 

.827 

.779 

.057 

76 

MD 

.820 

.751 

.058 

87 

.824 

.777 

.057 

76 

LG 

.790 

.724 

.054 

96 

. 804 

.759 

.054 

84 

MML 

.792 

.726 

.054 

95 

.796 

.752 

.053 

87 

RE 

.808 

.736 

.055 

92 

.803 

.762 

.054 

86 

MN 

.827 

.758 

.057 

84 

.825 

.779 

.056 

76 
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TABLE  II.  10  (Continued) 


Case  10a:  k=5,  p. =0.950,  n.=15,  i=l,2,...,5.  R =0.774 

FIRST  METHOD:  Fbar  = 3.7,  Nfail  = 14 


********gQ$ 

CL********** 

********qq$ 

Qjj*  ********* 

A.90 

m 

s 

ACL 

> 

• 

00 

o 

m 

s 

ACL 

ML 

.851 

.650 

.123 

76 

.799 

.693 

.116 

72 

MD 

.820 

.634 

.114 

88 

.789 

.680 

.112 

72 

LG 

.631 

.542 

.077 

97 

.701 

.626 

.091 

95 

MML 

.716 

.567 

.094 

97 

.699 

.618 

.094 

89 

RE 

.723 

.597 

.099 

95 

.727 

.647 

.098 

90 

MN 

.765 

.621 

.099 

89 

.765 

.667 

.100 

76 

Case 

10b: 

Same  as  10a 

SECOND  METHOD:  Fbar  = 3. 

8,  Ncorr  = : 

L145 

ML 

.831 

.643 

.121 

86 

.788 

.687 

.114 

72 

MD 

.802 

.628 

.112 

89 

.778 

.680 

.111 

76 

LG 

.619 

.538 

.076 

96 

.704 

.621 

.091 

93 

MML 

.700 

.565 

.091 

97 

.699 

.616 

.092 

89 

RE 

.716 

.595 

.097 

96 

.724 

.645 

.096 

90 

MN 

.801 

.633 

.110 

89 

.770 

.674 

.107 

79 
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TABLE  II. 11  (Continued) 


Case  11a:  Same  as  10a,  except  n^=30,  i=l,2,...5. 

FIRST  METHOD:  Fbar  * 7.6,  Nfail  = 0 

********90%  cl**********  ********gg%  cl********** 


A.90 

m 

s 

ACL 

> 

• 

00 

o 

m 

s 

ACL 

ML 

.794 

.681 

.084 

83 

.787 

.712 

.081 

76 

MD 

.782 

.673 

.081 

88 

.783 

.708 

.080 

76 

LG 

.729 

.633 

.071 

98 

.750 

.681 

.074 

88 

MML 

.737 

.636 

.073 

96 

.738 

.672 

.072 

89 

RE 

.742 

.652 

.075 

95 

.745 

.687 

.074 

89 

MN 

.776 

.677 

.078 

88 

.778 

.707 

.076 

77 

Case  lib:  Same  as  11a 

SECOND  METHOD:  Fbar  * 7.6,  Ncorr  = 551 


ML 

.791 

.680 

.084 

86 

.786 

.711 

.081 

/ 6 

MD 

.779 

.672 

.081 

88 

.781 

.707 

.079 

76 

LG 

.726 

.632 

.071 

98 

.748 

.680 

.073 

88 

MML 

.737 

.636 

.072 

96 

.738 

.672 

.072 

89 

RE 

.742 

.652 

.075 

95 

.745 

.687 

.074 

89 

MN 

.786 

.680 

.080 

88 

.781 

.708 

.078 

76 
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TABLE  11.12  (Continued) 


Case  12a:  Same  as  10a,  except  n^=50,  i*l,2,...,5. 

FIRST  METHOD:  Fbar  = 12.5,  Nfail  ■ 0 

********90%  CL*********  *********80%  CL********* 


A.90 

m 

s 

ACL 

o 

00 

• 

< 

m 

s 

ACL 

ML 

.785 

.703 

.061 

89 

.770 

.727 

.059 

81 

MD 

.778 

.698 

.060 

89 

.768 

.725 

.059 

81 

LG 

.750 

.676 

.056 

97 

.750 

.710 

.056 

89 

MML 

.750 

.676 

.056 

95 

.745 

.704 

.055 

89 

RE 

.754 

.686 

.056 

95 

.760 

.713 

.056 

86 

MN 

.782 

.706 

.058 

89 

.770 

.728 

.057 

81 

Case  12b:  Same  as  12a 

SECOND  METHOD:  Fbar  = 12.5,  Ncorr  = 187 


784 

.703 

.061 

89 

.770 

.727 

.059 

81 

778 

.689 

.059 

89 

.768 

.725 

.058 

81 

749 

.675 

.056 

97 

.750 

.709 

.056 

89 

750 

.676 

.056 

95 

.745 

.704 

.055 

89 

754 

.686 

.056 

95 

.760 

.713 

.056 

86 

785 

.706 

.059 

89 

.770 

.728 

.058 

81 

TABLE  11.13  (Continued) 


Case  13a:  k-5,  pi=0.947,  N^=15,  i=l,...,5.  Rg=0.762 

FIRST  METHOD:  Fbar  = 3.9,  Nfail  = 7 

********90%  CL**********  *********80%  CL********* 


A.  90 

m 

s 

ACL 

o 

00 

• 

< 

m 

s 

ACL 

ML 

.762 

.641 

.120 

77 

.799 

.684 

.114 

63 

MD 

.738 

.626 

.111 

89 

.789 

.677 

.110 

74 

LG 

.627 

.535 

.075 

98 

.701 

.619 

.089 

97 

MML 

.648 

.559 

.091 

98 

.699 

.610 

.092 

91 

RE 

.714 

.590 

.097 

96 

.729 

.640 

.096 

89 

MN 

.721 

.614 

.098 

91 

.765 

.660 

.098 

77 

Case  13b:  Same  as  13a 

SECOND  METHOD:  Fbar  = 4.0,  Ncorr  =1122 


ML 

.750 

.634 

.117 

87 

.788 

.677 

.111 

74 

MD 

.727 

.619 

.108 

89 

.778 

.671 

.107 

74 

LG 

.617 

.533 

.076 

95 

.704 

.614 

.089 

93 

MML 

.648 

.557 

.089 

98 

.699 

.608 

.089 

91 

RE 

.707 

.588 

.094 

98 

.726 

.638 

.094 

92 

MN 

.731 

.624 

.106 

79 

.770 

.665 

.104 

75 
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TABLE  II . 14  (Continued) 


Case  14a:  Same  as  13a,  except  n^=30,  i>l,2,...,5. 

FIRST  METHOD:  Fbar  « 8.0,  Nfail  - 0 

********90%  CL**********  **»*****80%  CL********** 


A.  90 

m 

s 

ACL 

A.  80 

m 

s 

ACL 

ML 

.794 

.668 

.085 

88 

.785 

.700 

.081 

78 

MD 

.782 

.661 

.081 

88 

.780 

.697 

.080 

77 

LG 

.729 

.623 

.072 

96 

.748 

.671 

.075 

88 

MML 

.737 

.626 

.073 

96 

.738 

.661 

.072 

89 

RE 

.743 

.641 

.075 

94 

.743 

.677 

.074 

86 

MN 

.776 

.666 

.078 

88 

.771 

.696 

.077 

79 

Case 

14b: 

Same  as 

14a 

SECOND  METHOD:  Fbar  » 8. 

0,  Ncorr  » 493 

ML 

.791 

.668 

.084 

88 

.782 

.699 

.081 

78 

MD 

.779 

.660 

.081 

88 

.778 

.696 

.080 

77 

LG 

.726 

.622 

.072 

98 

.745 

.670 

.074 

88 

MML 

.737 

.626 

.073 

96 

.738 

.661 

.072 

89 

RE 

.743 

.641 

.075 

94 

.743 

.677 

.074 

86 

MN 

.786 

.668 

.080 

88 

.778 

.697 

.078 

78 

i 

I 
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TABLE  11.15  (Continued) 


Case  15a:  Same  as  13a,  except  n^=50,  i*l,2,...,5. 

FIRST  METHOD:  Fbar  = 13.2,  Nfail  * 0 

********90%  CL**********  *********go%  CL********** 


A.90 

m 

s 

ACL 

o 

00 

• 

< 

m 

s 

ACL 

ML 

.766 

.690 

.060 

85 

.768 

.715 

.059 

77 

MD 

.760 

.685 

.059 

92 

.766 

.713 

.058 

77 

LG 

.733 

.664 

.055 

97 

.749 

.697 

.056 

85 

MML 

.734 

.665 

.055 

97 

.741 

.692 

.054 

86 

RE 

.743 

.674 

.056 

95 

.745 

.701 

.055 

85 

MN 

.767 

.693 

.058 

86 

.768 

.715 

.057 

77 

Case  15b:  Same  as  15a 

SECOND  METHOD:  Fbar  *13.2,  Ncorr  « 169 


ML 

.766 

.690 

.060 

85 

.768 

.715 

.059 

77 

MD 

.760 

.685 

.059 

92 

.766 

.713 

.058 

77 

LG 

.733 

.663 

.055 

97 

.748 

.697 

.056 

85 

MML 

.734 

.665 

.055 

97 

.741 

.692 

.054 

86 

RE 

.743 

.674 

.056 

95 

.745 

.701 

.055 

85 

MN 

.767 

.694 

.058 

85 

.768 

.716 

.057 

77 
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TABLE  11.16  (Continued) 


Case  16a:  k=15,  p^=0.995,  n^=5,  i=l,2,...,5.  Rg=0.928 

FIRST  METHOD:  Fbar  = 1.3,  Nfail  = 157 


********90% 

it  ******* 

********80% 

CL********** 

A.90 

m 

s 

ACL 

> 

• 

00 

o 

m 

S 

ACL 

ML 

.907 

.831 

.078 

68 

.928 

.860 

.071 

68 

MD 

.872 

.801 

.072 

68 

.912 

.847 

.068 

68 

LG 

.708 

.632 

.057 

68 

.818 

.748 

.058 

68 

MML 

.778 

.707 

.064 

68 

.824 

.756 

.063 

68 

RE 

.839 

.745 

.072 

100 

.850 

.792 

.069 

100 

MN 

.810 

.757 

.053 

68 

.855 

.803 

.053 

68 

Case  16b:  Same  as  16a 

SECOND  METHOD:  Fbar  * 1 

ML  .939  .827  .094 

.4,  Ncorr  = i 

68  .954 

6952 

.856 

.084 

68 

MD 

.901 

.798 

.087 

68 

.938 

.843 

.081 

68 

LG 

.646 

.620 

.030 

68 

.797 

.744 

.053 

68 

MML 

.778 

.701 

.064 

68 

.824 

.750 

.063 

68 

RE 

.839 

.740 

.072 

100 

.850 

.787 

.069 

100 

MN 

.895 

.804 

.079 

68 

.920 

.836 

.074 

68 

46 
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TABLE  11.17  (Continued) 


1 


l 

- 


Case 

17a: 

Same  as 

16a , 

except 

n.=30. 

i=l,2. 

• • • § 15 

FIRST 

1 METHOD:  Fbar  = 2.3,  Nfail  - 37 

********90%  CL***********  ***** 

***80% 

cl  ********** 

A.90 

m 

s 

ACL 

> 

• 

00 

o 

m 

s 

ACL 

ML 

.925 

.867 

.058 

68 

.939 

.887 

.053 

68 

MD 

.907 

.852 

.055 

92 

.932 

.880 

.052 

68 

LG 

.796 

.764 

.035 

92 

.861 

.830 

.042 

92 

MML 

.842 

.796 

.049 

92 

.871 

.828 

.047 

92 

RE 

.880 

.816 

.052 

100 

.889 

.846 

.050 

97 

MN 

.877 

.840 

.043 

92 

.904 

.865 

.042 

92 

Case  17b:  Same  as  17a 

SECOND  METHOD:  Fbar  = 2 

ML  .918  .865  .061 

:.4,  Ncorr  =*  6426 

76  .933  .884 

.055 

68 

MD 

.901 

.850 

.058 

81 

.926 

.878 

.054 

87 

LG 

.816 

.777 

.039 

92 

.877 

.835 

.046 

92 

MML 

.837 

.795 

.048 

92 

.867 

.826 

.046 

92 

RE 

.876 

.815 

.051 

100 

.886 

.845 

.049 

97 

MN 

.906 

.857 

.054 

76 

.923 

.878 

.051 

76 

TABLE  11.19  (Continued) 


Case  19a:  k*10,  p. -0.990,  n.«15,  i»12,...,10.  R -0.904 

* X s 

FIRST  METHOD:  Fbar  » 1.6,  Nfail  = 110 

********90%  cl**********  ********go%  CL*********** 


A.  90 

m 

s 

ACL 

o 

CO 

• 

< 

m 

s 

ACL 

ML 

.907 

.801 

.090 

78 

.928 

.833 

.082 

78 

MD 

.872 

.774 

.083 

78 

.912 

.821 

.079 

78 

LG 

.708 

.618 

.056 

78 

.818 

.729 

.063 

78 

MML 

.778 

.684 

.072 

78 

.824 

.733 

.070 

78 

RE 

.820 

.719 

.079 

100 

.833 

.766 

.076 

100 

MN 

.810 

.736 

.062 

78 

.855 

.782 

.062 

78 

Case  19b:  Same  as  19a 

SECOND  METHOD:  Fbar  - 1.8,  Ncorr  = 4289 


ML 

.939 

.794 

.102 

78 

.954 

.826 

.091 

72 

MD 

.901 

.768 

.094 

78 

.938 

.815 

.088 

78 

LG 

.648 

.611 

.037 

88 

.797 

.725 

.061 

78 

MML 

.778 

.678 

.071 

78 

.824 

.728 

.069 

78 

RE 

.820 

.714 

.078 

100 

.833 

.762 

.075 

100 

MN 

.892 

.773 

.087 

78 

.919 

.808 

.082 

78 

TABLE  11.20  (Continued) 


Case  20a:  Same  as  19a,  except  N^=30,  i=l,2,...,10 
FIRST  METHOD:  Fbar  = 3.0,  Nfail  = 23 

********90%  cl*********  ********3Q%  cl********** 


A.90 

m 

s 

ACL 

> 

• 

00 

o 

m 

s 

ACL 

ML 

.925 

.839 

.066 

81 

.897 

.861 

.061 

59 

MD 

.907 

.825 

.064 

81 

.890 

.855 

.060 

80 

LG 

.796 

.751 

.043 

95 

.849 

.811 

.050 

95 

MML 

.842 

.773 

.057 

95 

.834 

.805 

.054 

95 

RE 

.869 

.791 

.059 

98 

.868 

.822 

.056 

93 

MN 

.877 

.818 

.053 

95 

.877 

.844 

.052 

81 

Case  20b:  Same  as  20a 

SECOND  METHOD:  Fbar  *3.1,  Ncorr  =3703 


ML 

.919 

.837 

.068 

81 

.892 

.859 

.062 

78 

MD 

.901 

.823 

.065 

83 

.886 

.853 

.061 

81 

LG 

.816 

.758 

.048 

95 

.844 

.813 

.054 

95 

MML 

.837 

.772 

.056 

95 

.834 

.804 

.054 

95 

RE 

.865 

.790 

m 

.058 

98 

.868 

.821 

.055 

94 

MN 

.905 

.831 

.062 

81 

.885 

.853 

.058 

81 

50 


I 


TABLE  11.21  (Continued) 


Case 

1 

21a: 

Same  as 

19a, 

except 

ni=50 , 

• • • 1 10 

FIRST  METHOD:  Fbar 

= 5. 

0,  Nfail  = 5 

■ 

i 

********90%  CL********** 

********  30% 

C l** ******** 

A.90 

m 

s 

ACL 

A.80 

m 

s 

ACL 

ML 

.925 

.853 

.048 

80 

.913 

.870 

.044 

76 

MD 

.915 

.844 

.046 

89. 

.909 

. 867 

.044 

76 

LG 

.856 

.809 

.039 

99 

.880 

.843 

.040 

94 

: 

MML 

.876 

.812 

.043 

96 

.874 

.836 

.041 

89 

RE 

.879 

.823 

.044 

97 

.881 

.846 

.042 

91 

MN 

.910 

.847 

.042 

89 

.902 

.865 

.040 

80 

: 

Case 

21b: 

Same  as 

21a 

SECOND  METHOD:  Fbar  = 5 

.1,  Ncorr  = 3019 

ML 

.923 

.852 

.048 

86 

.911 

.869 

.045 

76 

MD 

.913 

.843 

.047 

89 

.907 

.866 

.044 

76 

LG 

' 

.869 

.809 

.040 

99 

.883 

.843 

.041 

89 

MML 

.876 

.812 

.043 

96 

.874 

.836 

.041 

89 

RE 

.879 

.823 

.044 

97 

.881 

.846 

.042 

91 

MN 

.918 

.851 

.045 

86 

.908 

.868 

.042 

76 

TABLE  11.22  (Continued) 


Case 

22a: 

k=15 , p 

^=0.990 

, i-1 

1 2 $ • • . / 

14 ' Pl5 

=0.995; 

ni=15 

1—1 ! 2 t • 

. . ,15. 

v°* 

864 

FIRST  METHOD:  Fbar  = 2.3 

, Nfail  = 62 

*********90% 

Cl*** ****** 

********30% 

CL********** 

A.90 

m 

s 

ACL 

A.80 

m 

S 

ACL 

ML 

.907 

.755 

.110 

62 

. 879 

.791 

.102 

62 

MD 

.872 

.732 

.102 

87 

.866 

.781 

.098 

62 

LG 

.708 

.589 

.061 

87 

.742 

.695 

.076 

87 

MML 

.778 

.648 

.085 

87 

.765 

.698 

.084 

87 

RE 

.796 

.683 

.093 

99 

.815 

.731 

.090 

94 

MN 

.810 

.704 

.080 

87 

.813 

.749 

.080 

87 

Case 

22b: 

Same  as 

22a 

SECOND  METHOD:  Fbar  = 2. 

4,  Ncorr  = 6499 

ML 

.939 

.744 

.116 

76 

.859 

.781 

.106 

76 

MD 

.901 

.722 

.107 

87 

.846 

.771 

.103 

78 

LG 

.646 

.590 

.055 

90 

.755 

.692 

.077 

87 

MML 

.778 

.643 

.083 

87 

.751 

.693 

.082 

87 

RE 

.796 

.678 

.091 

99 

.805 

.727 

.089 

94 

MN 

.895 

.732 

.102 

87 

.841 

.768 

.097 

74 

TABLE  11.23  (Continued) 


Case 

23a: 

Same  as 

22a, 

except 

nj_=30, 

i=l,2 

9 • • • f 15 

FIRST  METHOD:  Fbar 

=4. 

3,  Nfail 

= 6 

********90%  CL 

******** ** 

********QQ% 

CL********** 

A.90 

m 

s 

ACL 

o 

00 

• 

m 

S 

ACL 

ML 

.877 

.793 

.072 

78 

.897 

.818 

.067 

66 

MD 

.862 

.781 

.069 

92 

.890 

.813 

.066 

78 

LG 

.772 

.719 

.050 

99 

. 837 

.774 

.056 

92 

MML 

.803 

.733 

.061 

99 

.834 

.767 

.059 

92 

RE 

.833 

.751 

.064 

95 

.839 

.784 

.062 

90 

MN 

.853 

.780 

.060 

79 

.877 

.807 

.059 

78 

Case  23b:  Same  as  23a 

SECOND  METHOD:  Fbar  = 4.3,  Ncorr  = 5644 


ML 

.872 

.789 

.072 

78 

.892 

.815 

.067 

76 

MD 

.857 

.777 

.069 

92 

.886 

.810 

.066 

76 

LG 

.789 

.723 

.056 

99 

. 842 

.775 

.059 

92 

MML 

.803 

.733 

.061 

99 

.834 

.767 

.059 

92 

RE 

.833 

.751 

.064 

96 

.839 

.784 

.062 

90 

MN 

.865 

.789 

.067 

79 

.886 

. .812 

.064 

76 

TABLE  11.24  (Continued) 


Case 

24a: 

Same  as 

22a, 

except 

n.=50, 

i*l/2 , 

. . . ,15 

FIRST  METHOD:  Fbar  = 7.3,  Nfail  = 0 

*********90%  CL*********  ********80% 

CL********** 

A.90 

m 

s 

ACL 

> 

• 

00 

o 

m 

s 

ACL 

ML 

.875 

.805 

.053 

85 

.869 

.825 

.050 

75 

MD 

.866 

.797 

.052 

85 

.865 

.822 

.050 

76 

LG 

.830 

.767 

.046 

99 

.843 

.801 

.047 

90 

MML 

.832 

.769 

.048 

94 

.834 

.794 

.046 

85 

RE 

.848 

.780 

.049 

94 

.843 

.805 

.047 

88 

MN 

.869 

.804 

.049 

85 

.865 

.823 

.047 

77 

Case  24b:  Same  as  24a 

SECOND  METHOD:  Fbar  = 1 

ML  .873  .804  .053 

’.3,  Ncorr  = 4528 

85  .868  .824 

.050 

75 

MD 

.864 

.796 

.052 

92 

.864 

.821 

.050 

81 

LG 

.829 

.767 

.047 

98 

.842 

.801 

.047 

92 

MML 

.832 

.769 

.048 

98 

.834 

.794 

.046 

93 

RE 

.848 

.780 

.049 

94 

.843 

.804 

.047 

88 

MN 

.873 

.807 

.050 

85 

.867 

.825 

.048 

75 

TABLE  11.25  (Continued) 


Case  25a:  k=15,  p^=0.990,  N^=15,  i*l,2, . . . ,15.  Rg=0.860 

FIRST  METHOD:  Fbar  = 2.3,  Nfail  = 46 

********90%  CL**********  *********80%  CL********** 


A.90 

m 

s 

ACL 

> 

• 

00 

o 

m 

s 

ACL 

ML 

.851 

.751 

.103 

65 

.879 

.787 

.095 

65 

MD 

.820 

.728 

.095 

91 

.866 

.111 

.091 

65 

LG 

.681 

.587 

.054 

91 

.742 

.692 

.069 

91 

MML 

.716 

.645 

.078 

91 

.765 

.695 

.077 

91 

RE 

.783 

.678 

.084 

99 

.804 

.727 

.082 

95 

MN 

.765 

.702 

.073 

91 

.813 

.747 

.074 

91 

Case  25b:  Same  as  25a 

SECOND  METHOD:  Fbar  = 2.5,  Ncorr  = 6460 


TABLE  11.26  (Continued) 


Case  26a:  Same  as  25a,  except  n^=30,  i=l,2,...,15 
FIRST  METHOD:  Fbar  = 4.5,  Nfail  = 6 

********90%  CL*********  ********80%  CL********** 


A.90 

m 

s 

ACL 

o 

00 

• 

< 

m 

s 

ACL 

ML 

.877 

.785 

.077 

81 

.858 

.810 

.072 

67 

MD 

.862 

.773 

.074 

82 

.853 

.805 

.070 

%80 

LG 

.772 

.713 

.054 

99 

.815 

.767 

.061 

93 

MML 

.80  3 

.726 

.065 

93 

.805 

.760 

.063 

93 

RE 

.832 

.744 

.067 

96 

.835 

.111 

.065 

89 

MN 

.853 

.773 

.065 

82 

.846 

.800 

.064 

81 

Case  26b:  Same  as  25a 

SECOND  METHOD:  Fbar  = 4.6,  Ncorr  = 5591 


ML 

.872 

.781 

.077 

81 

.855 

.807 

.072 

78 

MD 

.857 

.770 

.073 

87 

.849 

.802 

.070 

79 

LG 

.789 

.716 

.060 

99 

.811 

.768 

.063 

93 

MML 

.803 

.726 

.065 

99 

.805 

.760 

.063 

93 

RE 

.832 

.744 

.067 

96 

.835 

.777 

.065 

89 

MN 

.865 

.781 

.071 

81 

.851 

.805 

.068 

78 

56 


I 


TABLE  11.27  (Continued) 


Case  27a:  Same  as  25a,  except  n^=50,  i-l,2,...,15 
FIRST  METHOD:  Fbar  = 7.4,  Nfail  = 1 

********90%  CL*********  ********80%  CL********** 


> 

• 

VO 

o 

m 

S 

ACL 

o 

00 

• 

< 

m 

S 

ACL 

ML 

.875 

.802 

.057 

80 

.869 

.822 

.054 

74 

MD 

.866 

.795 

.055 

86 

.865 

.819 

.053 

74 

LG 

.823 

.765 

.049 

98 

.844 

.799 

.050 

86 

MML 

.832 

.767 

.051 

94 

.834 

.792 

.050 

86 

RE 

.841 

.778 

.052 

95 

.846 

.802 

.050 

89 

MN 

.869 

.802 

.052 

86 

.865 

.821 

.051 

74 

Case  27b:  Same  as  27a 

SECOND  METHOD:  Fbar  =7.5,  Ncorr  = 4577 


ML 

.873 

.801 

.056 

86 

.868 

.821 

.054 

74 

MD 

.864 

.794 

.055 

86 

.864 

.818 

.053 

74 

LG 

.829 

.765 

.050 

98 

.842 

.799 

.050 

86 

MML 

.832 

.767 

.051 

94 

.834 

.792 

.050 

86 

RE 

.841 

.777 

.052 

95 

.846 

.802 

.050 

89 

MN 

.873 

.805 

.054 

86 

.867 

.822 

.052 

74 

57 


~ * k,V 


TABLE  11.28  (Continued) 


Case 

28a: 

k-15, 

Pi-0.995 

, i-1 

,2,. ..,14 

' p15 

-0.850, 

ni*15 

i-1,2, . . . ,15.  R -0.792 

S 

FIRST  METHOD:  Fbar  - 3.2,  Nfail  » 16 

******** 90%  ££,********«  *****ft**gQ% 

QJj*  ********* 

A.90 

m 

s 

ACL 

A.80 

m 

s 

ACL 

ML 

.851 

.678 

.124 

72 

.799 

.719 

.116 

63 

MD 

.820 

.659 

.115 

82 

.789 

.711 

.112 

82 

LG 

.681 

.563 

.081 

97 

.719 

.648 

.093 

87 

MML 

.716 

.585 

.098 

97 

.699 

.636 

.097 

82 

RE 

.753 

.618 

.10  3 

97 

.757 

.668 

.102 

88 

MN 

.765 

.635 

.100 

97 

.765 

.683 

.100 

82 

Case  28b:  Same  , 

SECOND  METHOD:  ] 

ML  .828  .664 

as  28a 

Fbar  - 3. 

.123 

4,  Ncorr  - 6548 

82  .783  .707 

.116 

76 

MD 

.799 

.647 

.114 

82 

.774 

.699 

.112 

78 

LG 

.641 

.548 

.072 

98 

.698 

.635 

.088 

97 

MML 

.700 

.578 

.096 

97 

.699 

.629 

.096 

97 

RE 

.739 

.611 

.102 

97 

.750 

.662 

.101 

90 

MN 

.808 

.663 

.110 

82 

.773 

.701 

.107 

76 

> 


TABLE  11.29  (Continued) 

Case  29a:  Same  as  28a,  except  n^=30,  i=l,2,...,15 
FIRST  METHOD:  Fbar  - 6.6,  Nfail  = 1 

********90%  CL**********  ********80%  CL********** 


A.  90 

m 

s 

ACL 

o 

00 

• 

< 

m 

s 

ACL 

ML 

.796 

.702 

.081 

83 

.790 

.733 

.078 

80 

MD 

.783 

.693 

.078 

91 

.785 

.729 

.076 

80 

LG 

.731 

.651 

.067 

98 

.754 

.701 

.069 

91 

MML 

.737 

.653 

.071 

98 

.745 

.689 

.070 

92 

RE 

.761 

.670 

.073 

96 

.768 

.705 

.071 

89 

MN 

.784 

.694 

.073 

92 

.785 

.724 

.072 

80 

Case 

29b: 

Same  as 

29a 

SECOND  METHOD:  Fbar  = 6 

.7,  1 

Ncorr  = 

5985 

ML 

.792 

.699 

.081 

90 

.787 

.730 

.077 

80 

MD 

.780 

.690 

.078 

91 

.783 

.726 

.076 

80 

LG 

.726 

.648 

.066 

99 

.750 

.698 

.069 

92 

MML 

.737 

.652 

.071 

98 

.745 

.688 

.070 

92 

RE 

.761 

.669 

.072 

96 

.766 

.705 

.071 

89 

MN 

.792 

.706 

.076 

89 

.786 

.732 

.074 

80 

1 


TABLE  11.30  (Continued) 


Case 

30a: 

Same  as 

28a, 

except 

n.=50. 

1=1,2, 

. . . ,15 

FIRST 

1 METHOD:  Fbar  = 10.9,  Nfail  = 0 

********90%  CL**********  *********go% 

CL********* 

A.90 

m 

s 

ACL 

> 

• 

00 

o 

m 

s 

ACL 

ML 

.803 

.723 

.061 

87 

.800 

.747 

.059 

78 

MD 

.796 

.718 

.059 

88 

.797 

.745 

.058 

78 

LG 

.766 

.695 

.054 

98 

.777 

.729 

.054 

87 

MML 

.766 

.693 

.056 

97 

.768 

.721 

.055 

89 

RE 

.774 

.703 

.057 

95 

.775 

.731 

.056 

87 

MN 

.800 

.724 

.057 

88 

.793 

.746 

.056 

79 

Case  30b:  Same  as  30a 

SECOND  METHOD:  Fbar  = 11. 0,  Ncorr  = 

ML  .802  .722  .061  87  .798 

5499 

.746 

.058 

78 

MD 

.795 

.717 

.059 

89 

.796 

.744 

.058 

78 

LG 

.764 

.694 

.054 

98 

.778 

.728 

.054 

87 

MML 

.766 

.693 

.056 

97 

.768 

.720 

.055 

89 

RE 

.773 

.703 

.057 

95 

.775 

.730 

.056 

87 

MN 

.807 

.731 

.058 

87 

.801 

.751 

.056 

78 

60 


TABLE  11.31  (Continued) 

Case  31a:  k=13,  R =0.723 

S 

PiS  .995  .985  .979  .988  .982  .980  .967  .995  .970 
.995  .968  .980  .900 


n^=15,  i=l,2 , ... ,13. 

FIRST  METHOD:  Fbar  = 4.5,  Nfail  = 5 


********90%  CL 

********* 

********30% 

QJj**  ******** 

A.90 

m 

s 

ACL 

A.80 

m 

S 

ACL 

ML 

.762 

.602 

.123 

84 

.731 

.647 

.118 

68 

MD 

.738 

.589 

.114 

84 

.723 

.641 

.114 

84 

LG 

.598 

.510 

.081 

99 

.663 

.588 

.095 

89 

MML 

.648 

.531 

.094 

93 

.655 

.581 

.095 

84 

RE 

.685 

.560 

.098 

94 

.690 

.610 

.098 

88 

MN 

.721 

.583 

.102 

86 

.709 

.628 

.103 

84 

Case  31b:  Same  as  31a 

SECOND  METHOD:  Fbar  = 4. 

ML  .745  .591  .120 

7,  Ncorr  = ■ 

84  .718 

4779 

.637 

.115 

77 

MD 

.723 

.579 

.112 

91 

.711 

.631 

.111 

74 

LG 

.607 

.502 

.080 

99 

.648 

.579 

.093 

93 

MML 

.648 

.526 

.094 

93 

.655 

.577 

.095 

86 

RE 

.685 

.555 

.098 

95 

.689 

.605 

.099 

88 

MN 

.735 

.593 

.110 

84 

.711 

.633 

.108 

77 

61 


TABLE  11.32  (Continued) 


Case 

32a: 

Same  as 

31a, 

except 

ni=30 

, i«l,2. 

...  ,13 

FIRST  METHOD:  Fbar 

= 9. 

3,  Nfail  * 0 

********90%  CL********* 

********80% 

CL********** 

A.90 

m 

s 

ACL 

A.80 

m 

s 

ACL 

ML 

.728 

.631 

.080 

85 

.727 

.663 

.078 

75 

MD 

.718 

.625 

.077 

90 

.723 

.661 

.077 

82 

LG 

.675 

.590 

.069 

96 

.695 

.637 

.072 

88 

MML 

.678 

.594 

.069 

96 

.682 

.630 

.069 

86 

RE 

.701 

.609 

.071 

93 

.704 

.644 

.071 

86 

MN 

.725 

.633 

.075 

87 

.724 

.663 

.074 

79 

Case 

32b  : 

Same  as 

32a 

SECOND  METHOD:  Fbar  = 9 

.4,  Ncorr  = 

3673 

ML 

.725 

.629 

.080 

87 

.724 

.661 

.077 

78 

MD 

.716 

.623 

.077 

90 

.721 

.659 

.076 

84 

LG 

.672 

.588 

.069 

97 

.692 

.635 

.071 

90 

i 

MML 

.678 

.593 

.069 

96 

.682 

.629 

.069 

86 

RE 

.701 

.608 

. 073 

93 

.704 

.644 

.071 

86 

MN 

.729 

.637 

.076 

85 

.725 

.665 

.075 

76 

TABLE  11.33  (Continued) 


Case  33a:  Same  as  31a, except  n^=50,  i*l,2, . . . , 13 
FIRST  METHOD:  Fbar  » 15.6,  Nfail  * 0 


****it****gQ% 

££******** 

********30% 

CL********** 

A.  90 

m 

s 

ACL 

o 

00 

• 

< 

m 

s 

ACL 

ML 

.729 

.652 

.063 

84 

.733 

.678 

.061 

77 

MD 

.724 

.648 

.061 

87 

.731 

.676 

.061 

77 

LG 

.700 

.628 

.058 

95 

.714 

.662 

.058 

84 

MML 

.699 

.630 

.057 

93 

.707 

.658 

.057 

84 

RE 

.710 

.639 

.058 

93 

.716 

.667 

.057 

82 

MN 

.733 

.659 

.060 

84 

.734 

.681 

.059 

77 

Case 

33b: : 

Same  as 

33a 

SECOND  METHOD:  Fbar  *15 

.6, 

Ncorr  = 

2752 

ML 

.729 

.652 

.063 

84 

.732 

.677 

.061 

77 

MD 

.723 

.648 

.061 

90 

.730 

.676 

.061 

77 

LG 

.699 

.628 

.058 

95 

.713 

.661 

.058 

84 

MML 

.699 

.630 

.057 

93 

.707 

.658 

.057 

84 

RE 

.709 

.639 

.058 

93 

.716 

. 666 

.057 

82 

MN 

.735 

.660 

.061 

84 

.736 

.682 

.060 

77 

63 


TABLE  11.34  (Continued) 


f 


Case  34a:  k=13.  R =0.723 

s 


995 

.985 

.979 

.988 

.982 

.980 

.967 

.995 

970 

.995 

.968 

.980 

.900 

150 

90 

75 

100 

125 

18 

28 

125 

63 

125 

59 

5 

19 

first  METHOD:  Fbar  = 15.5,  Nfail  = 0 


CL  ********** 

CL ★*★★★★*★★** 

A.90 

m 

s 

ACL 

O 

00 

• 

< 

m 

s 

ACL 

ML 

.749 

.621 

.110 

81 

.740 

.656 

.101 

72 

MD 

- 

- 

- 

- 

- 

- 

- 

- 

LG 

.721 

.594 

.108 

91 

.724 

.641 

.099 

80 

MML 

.721 

.594 

.109 

90 

.716 

.631 

.102 

81 

RE 

.737 

.611 

.104 

86 

.733 

.647 

.097 

77 

MN 

.491 

.404 

.065 

100 

.549 

.484 

.070 

100 

•Case  34b’:  Same  as  34a 

SECOND  METHOD:  Fbar  =15 

ML  .723  .603  .102 

. 6 , Ncorr  = 

90  .720 

2425 

.640 

.095 

80 

MD 

.656 

.559 

.086 

99 

.699 

.624 

.089 

88 

LG 

.683 

.570 

.097 

97 

.700 

.624 

.091 

37 

MML 

.685 

.571 

.098 

97 

.687 

.611 

.094 

91 

RE 

.703 

.590 

.095 

94 

.708 

.630 

.090 

85 

MN 

.553 

.453 

.078 

100 

.594 

.522 

.080 

100 

TABLE  11.35  (Continued) 


Case 

35  a': 

TV 

II 

1— 1 
in 
• 

R =0.860 
s 

Pi: 

.990  .990 

.990 

.990  .990  .990 

.990 

.990 

.990  .990 

.990 

.990  .990  .990 

.990 

ni: 

250  40 

120 

15  130  65 

70 

75 

100  90 

60 

60  20  30 

40 

FIRST  METHOD:  Fbar  * 11.6 

, Nfail  = 0 

********90% 

CL********* 

********80%  CL********** 

A.  90 

m 

s 

ACL 

> 

OO 

o 

3 

s 

ACL 

ML 

MD 

LG 

.881 

.803 

.064 

80 

.874  .823 

.059 

72 

.861 

.776 

.067 

89 

.862  .809 

.060 

79 

MML 

.862 

.778 

.068 

89 

.856  .802 

.062 

82 

RE 

.870 

.788 

.066 

86 

.863  .812 

.060 

78 

MN 

.780 

.720 

.049 

100 

.803  .762 

.048 

99 

Case 

35tt:  : 

Same  , 

as  35a 

SECOND  METHOD : 

Fbar  = 11. 

6,  Ncorr  = 3881 

ML 

.876 

.800 

.063 

82 

.870  .820 

.058 

73 

MD 

.855 

.784 

.060 

91 

.865  .815 

.057 

76 

LG 

.854 

.772 

.065 

92 

.857  .806 

.058 

81 

MML 

.856 

.774 

.066 

91 

.853  .799 

.061 

83 

RE 

.864 

.785 

. 064 

88 

.859  .809 

.059 

80 

MN 

.811 

.743 

.056 

99 

.824  .778 

.053 

95 
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TABLE  11.36  (Concluded) 


Case 

36a: 

k=15 . 

R =0.792 
s 

Pi  •*  • 

995  .995 

.995 

.995  . 

995  .995 

.995  . 

995 

• 

995  .995 

.995 

.995  . 

995  .995 

.850 

ni : : 

20  20 

20 

20 

20  20 

20 

20 

20  20 

20 

20 

20  20 

150 

FIRST  METHOD::  Fbar  = 23. 

9,  Nf ail  = 

0 

******** 9Q% 

Ql^********* 

*********gQ% 

QT  ********* 

A.90 

m 

s 

ACL 

A.80 

m 

s 

ACL 

ML 

.824 

.732 

.072 

78 

.808 

.754 

.065 

71 

MD 

- 

- 

- 

- 

- 

- 

- 

- 

LG 

.814 

.714 

.075 

81 

.803 

.744 

.067 

75 

MML 

.812 

.713 

.075 

81 

.798 

.737 

.068 

77 

RE 

.816 

.721 

.072 

81 

.804 

.745 

.066 

75 

MN 

.733 

.677 

.048 

99 

.751 

.713 

.049 

96 

Case 

36b:. 

Same 

as  36a 

SECOND  METHOD : 

Fbar  = 24 

.3,  Ncorr  = 

6296 

ML 

.802 

.711 

.069 

86 

.791 

.735 

.064 

80 

MD 

.788 

.697 

.067 

91 

.785 

.730 

.063 

81 

LG 

.790 

.692 

.071 

88 

.781 

.724 

.064 

84 

MML 

.789 

.691 

.071 

88 

.776 

.718 

.066 

84 

RE 

.794 

.701 

.069 

88 

.785 

.728 

.064 

83 

MN 

.749 

.676 

.058 

97 

.759 

.709 

.057 

93 

66 


APPENDIX 


THE  COMPUTER  PROGRAM 

The  computer  program  consists  of  two  main  programs  and 
five  subroutines . The  two  main  programs  are  basically  the 
same.  The  first  main  program  is  used  with  the  FIRST  METHOD 
of  introducing  partial  component  failures  and  the  second  main 
program  is  used  with  the  SECOND  METHOD  of  introducing  partial 
component  failures.  The  main  program  reads  the  input  (case) 
k,  n^,  pi,  i=l,2,...,k.  With  this  input  the  main  program 
constructs  binomial  data  and  sends  this  data  to  subroutines 
MADSKY  (Madansky) , RMMLI  (Easterling  and  Randomized  Easterling) , 
AMLMAN  (Max.  I'hood  and  Mann)  and  WOODBG  (Log-Gamma);  these 
subroutines  compute  lower  confidence  limits . For  each  case 
these  computations  are  repeated  500  times.  Then  subroutine 
CMPARE  computes  the  90-th  and  80-th  percentile,  mean,  std. 
dev.  and  actual  confidence  level. 

The  programs  are  self-explanatory  since  comments  are 
inserted  throughout  the  programs . Note  that  subroutine  RANDOM 
and  PXSORT  are  from  Ref.  10,  subroutine  MDBETI  from  Ref.  3. 

The  percentile  points  c,f  the  chi-square  and  standard  normal 
distributions  are  from  Ref.  8. 

The  successful  root  finding  technique  in  the  MD  procedure 
which  was  developed  by  Lisowsky  (Ref.  4)  is  implemented  in 


subroutine  MADSKY. 
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